Josephson effect in mesoscopic graphene strips with finite width 
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We study Josephson effect in a ballistic graphene strip of length L smaller than the superconduct- 
ing coherence length and arbitrary width W . We find that the dependence of the critical supercurrent 
Ic on W is drastically different for different types of the edges. For smooth and armchair edges 
at low concentration of the carriers 7c decreases monotonically with decreasing W/ L and tends to 
a constant minimum for a narrow strip W/L < 1. The minimum supercurrent is zero for smooth 
edges but has a finite value eAo/h for the armchair edges. At higher concentration of the carriers, 
in addition to this overall monotonic variation, the critical current undergoes a series of peaks with 
varying W. On the other hand in a strip with zigzag edges the supercurrent is half-integer quantized 
to (n -I- l/2)4eAQ/?i, showing a step-wise variation with W . 

PACS numbers: 74.45.-(-c, 74.50.-fr, 73.63.-b, 74.78.Na 



Over the past years fabrication of the nanoclcctronic 
structures has provided the possibility of detecting the ef- 
fects of the electronic transport through a few or even sin- 
gle quantum states. The generic effect is the quantization 
of the conductance of a quantum point contact (QPC) 
to 2e^ /fi times the number of the transparent quantum 
channels^'i^. In a superconducting quantum point contact 
(SQPC) the analogous behavior occurs for the Josephson 
current*^. Furusaki et alA have predicted that the criti- 
cal (maximum) supercurrent through a two-dimensional 
electron gas (2DEG) junction shows a step- like variation 
as a function of the width of the junction depending on 
the Fermi wave length of the electrons. The quantization 
of the supercurrent in units of eAo/h (Aq is the super- 
conducting gap in the electrodes) was confirmed experi- 
mentally in semiconducting structures^. 

Recently experimental realization of graphene, the two 
dimensional carbon atoms arranged in a honeycomb lat- 
tice, has introduced a new type of mesoscopic mate- 
rial with unique properties^iii^. Most of the peculiar- 
ities come from the electronic structure of graphene 
which is fundamentally different from that of a metal 
or a semiconductor. Graphene has a gapless semi- 
metallic band structure with a linear dispersion rela- 
tion of the low-lying excitations. This makes the elec- 
trons in graphene to behave identical to two dimensional 
massless Dirac Fermions^iiSiiii. Already several quantum 
transport phenomena including the integer quantum Hall 
effeclji'2'i^ , the conductance quantization^^ and quantum 
shot noisei^ have been revisited in graphene and found 
to have anomalous features due to the relativistic-like 
dynamic of the electrons. 

Very interestingly graphene as a weak link between 
two superconductors can serve as a relativistic Joseph- 
son contacti^. Together with the fact that the graphene 
samples are realized with different type of the edges sug- 
gest a new class of Josephson devices with novel proper- 
ties. In particular the quantization of the supercurrent 
through a graphene strip with finite width is expected to 
depend on the type of the edges. The aim of this paper 
is to study graphene Josephson junctions with different 
edges. 



Josephson effect in a wide graphene sheet contact- 
ing two superconductors was studied by Titov and 
Beenakker^''^ just recently. Using a Dirac-Bogoliubov- 
de Gennes (DBdG) formalism^^, they found that in a 
ballistic graphene a Josephson current can flow even in 
the limit of zero concentration of the carriers i. e. at the 
Dirac point. The dependence of the critical current on 
the junction length as well as the current-phase relation 
were found to be the same as in a Josephson junction 
with an ordinary disordered normal meta l^^d^ . 

In this paper using the formalism of Ref[l^ we study 
the Josephson effect in a graphene strip with length 
L smaller than the superconducting coherence length 
^ = hv / Aq and an arbitrary width W for smooth, arm- 
chair and zigzag edges. We find that in contrast to an 
ordinary SQPC^i^ the supercurrent in smooth and arm- 
chair strips with a low concentration of the carriers is not 
quantized but rather shows a monotonic decrease with 
lowering W . For a narrow strip W L the supercurrent 
reduces to a constant minimum which is eAo/h for arm- 
chair edges and vanishingly small for smooth edges. Far 
from the Dirac point at a high concentration of the carri- 
ers the behavior of the supercurrent versus W consists of 
a similar monotonic variation and a series of peaks with 
distances inversely proportional to the chemical potential 
fi. We explain absence of the supercurrent quantization 
as a result of the significant contribution of the evanes- 
cent modes in the supercurrent which is a unique prop- 
erty of graphene monolayers. For a zigzag strip the situ- 
ation is drastically different because of the valley filtering 
characteristioi^ of the wave functions. We find that this 
type of the edges supports a half-integer quantization of 
the supercurrent to {n + l/2)4eAo/?i. 

The geometry of the studied Joseohson junction is 
shown schematically in Fig. [T] A ballistic graphene strip 
(N) connects two wide superconducting regions (S) which 
are produced by depositing of superconducting electrods 
on top of the graphene sheets. We assume the interfaces 
between the superconductors and the graphene strip to 
be free of defects and impurities. The superconduct- 
ing parts are assumed to be heavily doped such that 
the Fermi wave length Xp s inside them is very smaller 
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FIG. 1: Schematic of the studied junction structures. Right 
(left) junction has armchair (zigzag) edges. 



than the superconducting coherence length ^ and also 
Fermi wave length of the normal graphene strip Xfn- 
By the first condition mean field theory of superconduc- 
tivity will be justified and by second we can neglect the 
spatial variation of A in the superconductors close to 
the normal-superconductor (NS) interfaces. Thus the su- 
perconducting order parameter has the constant values 
A = Ao exp(±«</)/2) in the left and right superconductor, 
respectively, and vanishes identically in N. 

The superconducting correlations between electron- 
hole excitations are described by DBdG equation for the 
electron and hole wave functions u and 
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Here e is the excitation energy and the single particle 
Dirac Hamiltonian H = ~ihv{dx ctiTq + dya2T3), the 
time-reversal operator Q — cs ti , and the superconduc- 
tivity pair potential A = A ctq ''"i : operate on the wave 
functions u and v. The wave functions are 4-component 
spinors {ipi, ■02, i^'i-, ^^2) in which the indices show the two 
sublattices of the honeycomb lattice and primed (un- 
primed) V' indicates K' (K) Dirac point in inverse hon- 
eycomb lattice. CTi and (i — 1,2,3), are Pauli matrices 
in the sublattice and valley spaces respectively with the 
corresponding unit matrices ctq and tq . Eq. [T] describes 
the superconducting correlations between massless Dirac 
electron and hole excitations with opposite spins and dif- 
ferent valleys. 

The phase difference (p between the order parameters 
drives a supercurrent through the graphene strip which 
constitutes a weak link between two superconductors. 
This Josephson supercurrent is carried by the so called 
Andreev states which are formed in N region due to the 
successive conversion of the electron-hole excitations to 
each other (Andreev reflection) at the NS boundaries. 
For a short junction of L <C ^ and at zero tempera- 
ture the Andreev bound (discrete) states with energies 
e| < Ao have the main contribution to the Josephson 
supercurrent^. In this case we can neglect the contri- 
bution from the continuous states above the supercon- 
ducting gap. We obtain the energies e{(j)) of Andreev 
bound states by solving the DBdG equations with ap- 



propriate boundary conditions which are described be- 
low. The Josephson current is then obtained from the 
formula 



4e 



(2) 



where the factor 4 accounts for the spin and valley de- 
generacies. 

Inside N the solutions of DBdG equations are the elec- 
tron and hole-like wave functions which are classified by 
the 2-dimensional wave vector k = {k,q) with the energy- 
momentum relation e — hv\k.\. For a finite width W the 
transverse momentum is quantized (qn) by imposing the 
boundary conditions at the edges. The transverse bound- 
ary conditions are different for different edge types which 
a graphene strip can have. 

The solutions inside the S are rather mixed electron- 
hole excitations, called Dirac-Bogoliubov quasiparticles. 
In general to match the solutions inside different regions 
one needs the scattering matrices at the interfaces. But 
assuming ideal NS boundaries electrone-hole scattering 
can be described by a longitudinal boundary condition 
between electron and hole wave functions which readi^. 



/3 — arccos(£/Ao) 
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where tp is the phase of A in S and n is the unit vector 
perpendicular to the interface pointing from N to S. 

For each type of the edges by imposing the longitudinal 
boundary conditions at the interfaces we obtain the An- 
dreev energies en{4') of the transverse mode n. We will 
consider the regime Aq <^ /i where the retro-Andreev 
reflection dominates^^. In this limit the longitudinal mo- 
mentums of electron and hole fc„(±e) can be approxi- 
mated with fc„(0) . 

Let us start with analyzing a graphene strip with 
smooth edges which corresponds to an infinite mass con- 
finement. For this case the transverse boundary condi- 
tions do not mix the valleys and have the form V'i(O) = 
-V'2(0) and i^iiW) = ip2{W)^. For the electron-like An- 
dreev bound state the two-component wave function is a 
composition of the right and left going waves u{k) and 
u{—k) with u{k) = e''^^[sin(g?/ — a/2), sin(5y -fa/2)], and 
a quantized transverse momentum qn = (n + ^)n/W, 
[n = 0,1,2,...). The angle a = arctan(g/fc) indicates 
the propagation direction. Correspondingly the hole-like 
wave functions are written as a composition of v{k) and 
v{—k) which are the time-reversed of the electron wave 
functions. Imposing the electron- hole conversion condi- 
tion (Eq. [3]) at the NS boundaries {x — 0, L) we obtain 
the energy of the Andreev state as 



,(0) = Aoyi 



<„sin2(0/2), (4) 
~^ is the transmission 



where t„ = (g„/fc„)2 siv?{knL)] , 

coefficient of the mode n through the junction 
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For the armchair edges the two valleys are mixed by 
the boundary conditions ?Ai = and ip2 = V'2 
two edges^i. The electron wave function is a mixture of 
the two valleys with the form 

u = e'''^[e">y{e-"^, e"), (-l)"e-'««(e"*", e*")], (5) 

where the transverse momentum quantized to g„ = 
mr/W with integer n. Note that in contrast to the 
smooth edges, the wave functions u(q) and u(—q) rep- 
resent distinct solutions in which the components of the 
two valleys are interchanged. As a result the lowest mode 
n = has not the two fold valley degeneracy of the higher 
n = 1, 2, ... modes. Like the case of smooth edges the elec- 
tron wave function of a given mode n is the composition 
of u(fc) and u{—k). We find that the Andreev energies 
have the same form as in the smooth strip (Eq. [4]). 

We now consider the case of zigzag edges. In contrast 
to the smooth and armchair edges, these have a property 
of coupling the longitudinal and transverse momentums 
k and q. The boundary conditions are = ^'-^ ^ Q dX 
y = and ?/)2 = = y = W^^ , which lead to 
the transcendental relation sm{qW)/qW ~ ±hv/ fxW for 
the transverse momentum q. This relation has a finite 
number of solutions depending on the value of ^W/fiv. 
For ^W/hv < 1 there is an imaginary solution q labeling 
an evanescent mode in the y-direction. In the interval 
1 < ^W/hv < 3tt/2, the transcendental relation has a 
single oscillatory solution. For larger jiW/hv the number 
of the solutions are increased by two whenever the width 
is increased by a half of Fermi wave length Xp n = hv/ fj,. 
We classify the wave functions to two groups according 
to ± sings in the transcendental relation. The electronic 
wave functions of the first group (-1- sign) have the forms 

u+ = e*''"' [sm{qy) , sin{qy -f- a) , 0, 0] , (6) 

u~ = e-^'=^[0,0,sin(q2/),sin(gy-a)], (7) 

which, respectively, describe a right-going wave on the 
valley K and a left-going wave in the other valley K'. 
The wave functions of the second group (— sign) is 
obtained from the above functions by the replacement 
k —k. We note that for each mode the zigzag strip 
operates as a valley filter for the waves with mixed valley 
components'^. Since the valley index and q are conserved 
upon a normal reflection, the filtering property prevents 
the electrons from such scatterings. Consequently the 
wave functions of the electron bound states in N will be 
a composition of and From the longitudi- 

nal condition (Eq. [3]) we obtain that the energy of the 
Andreev state does not depend on the transverse mode 
q and is given by e = Ag cos(0/2), which is the same as 
for an ordinary SQPC^. 

From Eq. [5] and using the above obtained results for 
the allowed transverse momentums q and the Andreev 
energies e we obtain the Josephson current for the three 




FIG. 2; The critical current vs. relative width of junction for 
two different values of jiL/hv. Solid (dashed) line is corre- 
sponding to junction with smooth (armchair) edges. 



edge types. The results of smooth and armchair edges 
are similar. In both cases there are infinite numbers of 
transverse modes. For a given ^L/Tiv there are oscillatory 
modes with real k as well as evanescent modes with imag- 
inary k. Unlike the ordinary QPCs the evanescent modes 
in the graphene strip may have a significant transmission 
depending on the value of ^L/Tiv. In particular at the 
Dirac point ^L/Tiv <C 1 there is only few propagating 
modes and the evanescent modes have the main contri- 
bution in the supercurrent. This is the unique property 
of the graphene Josephson junctions in which very small 
Fermi energies are accessible. 

Fig. [2] shows the critical current dependence on width 
for the smooth and armchair edges and for two values 
of fiL/Tiv. For a typical small value jiL/hv = 0.1 (Fig. 
[21d) Ic has a linear dependence on W/L at large W/L 
which is the result of a diffusive-like transport in the bal- 
listic graphenc^^. By decreasing W/L the transmission 
of the modes through N decreases and the critical cur- 
rent shows a monotonic decrease without any quantiza- 
tion. For a narrow strip W ^ L the evanescent modes 
transparencies is vanishingly small and Ic reduces to a 
constant minimum value. For the smooth edges no prop- 
agating mode can exist when W < {ir /2)'hv / ^. So the 
minimum supercurrent goes to zero. However in the case 
of the armchair edges a lowest nondegenerate mode with 
g = always can propagate irrespective of the width of 
the junction. This zero mode result in a bound state 
carrying a nonzero residual supercurrent Ic — eAo/?i for 
W <L (see Fig. W)- 

For a higher chemical potential, ^L/Tiv = 5 (Fig. [2^) 
an oscillatory variation is added to the monotonic behav- 
ior. In this regime both the propagating and evanescent 
modes contribute to the supercurrent. In addition to an 
overall increase with ^W/fiv the critical current under- 
goes peaks which become smoother by increasing ^W/Tiv. 
Each peak (jump) signals addition of a new propagating 
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lowing important anomalies (see Fig. [3]). In contrast to 
an ordinary SQPC''''*, the width of the first step (3/4) 
is bigger than that of the higher steps (1/2). The extra 
width is the contribution of the single y-dircction evanes- 
cent mode (imaginary q) for W < hv/fi. Also the height 
of the first step is 1/2 of the height of the higher steps 
which itself is four times bigger than eAo/h the super- 
current quanta in an ordinary SQPC. Therefore the su- 
percurrent through a zigzag graphene strip is half-integer 
quantized to (n -I- l/2)4eAo/7i. The effect resembles the 
conductance quantization in graphene strips with zigzag 
edgesi^ and also the half-integer Hall effeciji'^'i^ in mono- 
layer structures. 



FIG. 3: The critical current vs. relative width of junction 
with zigzag edges. 

mode in the transport. Note that due to a nonzero trans- 
verse momentum of the lowest mode (n — 0) in smooth- 
edge strip the oscillations are shifted by half a period 
with respect to those of the armchair-edge one. Overall 
absence of the supercurrent quantization in these junc- 
tions is a consequence of the existence of the evanescent 
modes with appreciable transmission through N. 

Finally let us analyze the Josephson supercurrent for 
the zigzag strip. As we argued above in this case the 
valley filtering nature of the wave functions in N pre- 
vents the normal reflections at NS-interfaces. There is 
no evanescent mode with imaginary k and the situation 
is similar to an ordinary SQPC. The critical supercurrent 
shows a step- wise variation with fiW/hv but with the fol- 



In conclusion we have investigated the Josephson effect 
in a short monolayer graphene strip connecting two su- 
perconducting electrodes. Within a Dirac-Bogoliubov-de 
Gennes formalism, we have found that the variation of 
the critical supercurrent Ic versus the width W is drasti- 
cally different for the graphene strips with different edges. 
In the smooth and armchair strips with low concentration 
of the carriers Ic decreases monotonically with decreas- 
ing W/L. For a narrow strip W ^ L, Ic takes a constant 
minimum value which is given by eAo/h for the arm- 
chair edges but vanishingly small for the smooth edges. 
Increasing the concentration of the carriers, this overall 
monotonic dependence of Ic acquires a series of peaks 
with distances inversely proportional to the chemical po- 
tential fj,. For the zigzag edges we have found a step- 
wise variation of Ic versus iiW/hv, implying a half-integer 
quantization of the supercurrent to (n4-l/2)4eAo/?i with 
integer n. 
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